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Within the framework of cosmological inflationary models, relic gravitational waves arise as a nat- 
ural consequence of the quantum nature of primordial space-time metric fluctuations and validity 
of general relativity theory. Therefore the detection of cosmological gravitational waves would 
be the strongest evidence in support of basic assumptions of inflation theory. Although tracing 
gravitational waves by polarization patterns on last scattering surface of the cosmic microwave 
background is only planned for forthcoming experiments, some general constraints on the tensor 
mode of metric perturbations (i.e. gravitational waves) can be established right now. 
We present the determination of the amplitude of relic gravitational waves power spectrum. In- 
direct best-fit technique was applied to compare observational data and theory predictions. As 
observations we have used data on large-scale structure of the Universe and anisotropy of cos- 
mic microwave background temperature. The conventional inflationary model with 11 parameters 
has been investigated, all of them evaluated jointly. This approach gave us a possibility to find 
parameters of power spectrum of gravitational waves along with statistical errors. 
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INTRODUCTION 



For the last decade due to the progress in observations the cosmology entered the new stage of its development. 
This stage was signalized by the new generation of experiments aimed to the measurements of anisotropies of 
cosmic microwave background (CMB). These were balloon-borne BOOMERanG, MAXIMA, Archeops, ground- 
based interferometers DASI, CBI, VSA. And probably the most important one is the successor of COBE space 
mission, Wilkinson Microwave Anisotropy Probe (WMAP), that published 1-year observation results in 2003. 
WMAP has carried out the measurements of CMB over the whole sky with unprecedented angular resolution 
and high sensitivity of detectors. The data and web-links for these experiments are available at the web-site 
of Legacy Archive for Microwave Background Data Analysis (LAMBDA) project p. CMB explorations were 
complemented with extensive measurements of expansion dynamics of the Universe by distances to Supernovae 
and large-scale structure surveys. 

Advances in quality of experimental data manifestly call for the model capable to explain the whole set of 
collected data. Now it's well understood, that the simplest cosmological models can not match the observations 
adequately as e.g. standard flat CDM model with scale-invariant power spectrum of density fluctuations. The 
observations advance the complication of model, the number of parameters increases as well as the number 
of phenomena encompassed by theory. Today elaborate inflationary models need about 11 parameters for the 
proper description of reality 

Cosmological parameters could be classified in a following way: 

• The parameters related to the background model of the homogeneous and isotropic Universe. Evolution 
of Universe in this model is determined by the amount of energy densities of different components in ratio 
to the critical one £li — pi/p cr - These are densities of baryon component f2;,, hot dark matter (massive 
neutrinos) cold dark matter Sl c dm, and density parameter for dark energy (i.e. cosmological constant 
or quintessence) fide- 



The global properties of Universe. According to the Friedmann equations sum of density parameters gives 
unity for the spatially flat Universe, or 1 — S\. for curved, so Qk is a curvature parameter. The Hubble 
constant Hq = lOO/i km/(s ■ Mpc) is regarded as a global parameter too. 
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• Parameters associated with the inflation. According to inflationary scenario the large-scale structure of 
Universe is assumed to be formed due to growth of primordial matter density perturbations because 
of the gravitational instability. The perturbations have adiabatic nature and originate from quantum 
fluctuations stretched to the cosmological scales during the stage of exponential expansion - inflation. 
Perturbations are usually described by their power spectra in most general form P s (k) = A s k Ua for scalar 
and Pt(k) — A t k nt for tensor mode of perturbations. 

• The list of parameters is concluded by some additional parameters, like late re-ionization parameter t c 
and biasing parameters to relate the distribution of density to the spatial distribution of real astrophysical 
objects, e.g. galaxies. 

This parameter set is required for building the predictions of model to be compared with the observations 
for the purpose of finding out how good the theoretical assumptions match the results of current observations. 
As far as these predictions depend on the parameters quite nonlinearly, the determination of the best-fit value 
for one parameter requests the determination of full set of parameters, by means of statistics. Actually that is 
the main task of this investigation, sometimes referred as a testing of cosmological models. 

In this paper we shall attract attention to the one parameter among many others, namely the primordial 
power spectrum of tensor mode of space-time metric perturbations (relic gravitational waves). We shall use the 
most common kind of cosmological model, and analyze some particular inflation models. Sometimes this model 
is designated as concordance model, or even standard model. The purpose of this model is to give plausible 
explanation for large-scale structure of the Universe and its global properties. 

THE NATURE OF RELIC GWs 

Indeed, gravitational waves are very simple conceptually but proved to be extremely elusive for detection. 
Gravitational waves ought to be emitted by any physical system with changes in quadrupole distribution of 
stress-energy density. This follows from equations of general relativity as a free solution (wave solution) for 
small linear space-time metric perturbations on the background metric. 

Obviously, that only astrophysical objects could produce gravitational waves of significant amplitudes. The 
number of probable astrophysical sources is under discussion for observational programmes (see for review 0). 
There is no direct successful detections by now of any kind of sources. For cosmological GWs in fact very small 
chance exist to be detected by any human-made antenna because of their super-long wavelengths. 

Cosmological GWs play the role of most ancient relic in our Universe. They have their origin in first 
moments in Universe evolution, presumably the times of inflation. If we assume general relativity theory to be 
valid in those times and quantum zero-point oscillations to exist, so we come to conclusion of GWs inescapable 
existence. Variable gravitational field in very early Universe parametrically amplifies quantum oscillations, 
making up stochastic gravitational background. The amplitude of these GWs is determined by energy scales 
at the moment of emission, so for inflation this amplitude determines the energy scale of processes that give 
rise for inflation. So the nature of relic gravitational waves is very fundamental for physics. In fact both scalar 
mode and tensor modes of space-time metric perturbations share common origin from single physical process, 
quantum oscillations. 

The subsequent evolution of metric linear perturbations is completely described by gauge- invariant formalism 
for cosmological perturbations. Relic GWs are represented by tensor mode of space-time metric perturbations 
on the isotropic and homogeneous expanding background. As far as GWs do not interact with rest of medium 
in expanding Universe, their amplitude should decrease with time (see [2] for exhaustive review). 

Gauge-invariant theory leads to equation for evolution of tensor perturbations of space-time metrics in 
vacuum or perfect fluid: 

H(T) +2 «H( T ) + (2K; + fc 2 )H( T >=0 (1) 
a 

where tP T ) is gauge-invariant amplitude of tensor perturbations for two polarizations, a - scale factor, K = 
(—1,0,1) - curvature index and k is the wavenumber. The dots denote the derivatives with respect to the 
conformal time r\. The solution of this equation is propagating damped wave. 



OUTLINE OF METHOD 

As it was stated above both tensor and scalar modes of space-time metric perturbations have the same 
origin. The perturbations of scalar mode connected to the density perturbations so they eventually lead to the 
formation of galaxies, clusters and voids, i.e. large-scale structure of the Universe. Since tensor perturbations 
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do not lead to the formation of large-scale structure, cosmological GWs can not be revealed by present state of 
observable structure. Both scalar and tensor perturbations produce the power in CMB angular power spectrum 
due to Sachs- Wolfe effect. Of course, because of its specific polarization properties relic GWs should generate 
particular polarization pattern of CMB anisotropics, but the detection of polarization patterns in CMB is matter 
of future. 




Figure 1. The scalar and tensor contributions to CMB angular power spectrum for particular model. 

Therefore we propose the indirect method. As far as CMB power spectrum consists of contributions from 
scalar and tensor modes, we can extract the last one if we know precisely from large-scale structure the scalar 
contribution. Also, relic GWs have extra-long wavelengths of particle horizon size, so the CMB power spectrum 
from them should manifest fast decrease with higher multipole moments (smaller scales). On the contrary, the 
density perturbations at larger multipoles produce a feature-rich CMB power spectrum, a series of acoustic 
peaks. Thus, the physical properties of tensor mode of cosmological perturbations of space-time metric can be 
estimated on the basis of observational data on: angular power spectrum of CMB combined and large-scale 
structure of Universe in wide range of scales - from galactic ones to the size of particle horizon, from 10 -3 
to 10 4 Mpc. This way we can constrain total level of intensity of gravitational waves. In multipoles range 
of i ~ 10 — 20 scalar and tensor contributions are comparable (Fig. So the precision of result is to be 
determined by precision of datapoints within this range. This method is basically a normalization of angular 
power spectrum. 

BUILDING PREDICTIONS 

Unfortunately, the cosmological models are quite complicated and the most of predictions could not be 
reduced to analytical functions with cosmological parameters as arguments. Therefore for the given parameter 
set one need to run numerical computations to make predictions. One way is the "brute-force" method, when 
the numerical packages like CMBfast 0] are utilized. Since the testing requires large number of models to 
be evaluated, so even quite fast codes need a huge computational resources to accomplish the task. Another 
demerit lies in so-called "black box" uncertainty of numerical calculations, when the physical processes are 
hidden inside computations. So we propose a semi-analytical approaches to make predictions of CMB power 
spectrum. 

For tensor contribution to CMB power spectrum within low multipoles range (2 < I < 20) we have devel- 
oped an analytical approximation able to reproduce computed spectrum with high accuracy, for wide range of 
parameters values. This approximation takes a form: 

t{i + l)Cj = r~~~T~~~~~p~T • Kt(Slde) x exp(-CK Slo) ■ e 2 + D(n t , fi„) • £) (2) 
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where the coefficients A, b, C and D represent polynomial functions of n t and fio = 1 — fife, K^dde) is 
amplification factor for CMB power spectrum caused by dark energy. Analytical formulae for them presented 
on our paper 

Semi-analytical approximation for scalar contribution to CMB power spectrum for the multipoles range 
2 < £ < 20 was developed in our previous paper |S] and used here. This approximation combine the speed 
of computations and clear unobscured physical meaning of processes that cause the anisotropies in CMB. The 
angular power spectrum of CMB at higher multipoles have a number of peaks separated with dips and manifest 
the damping of average amplitude to the highest I. This features are explained by acoustic oscillations in 
photon-baryon fluid set up by adiabatic perturbations on entering the sound horizon. Instead of calculations 
of power for each £ within a wide range, we propose to describe the whole shape of CMB power spectrum by 
positions and amplitudes of the peaks and dips. The insignificant loss of information will be a price for the 
computations speedup. The number of analytical approximations was developed for heights and positions of 
first three peaks and for position of first dip in CMB power spectrum in models under consideration, see [§]• 

OBSERVATIONAL DATA 

Of course, the results of testing of the model strongly depend on the data used. The next important 
requirement the data should meet is the statistical independence, i.e. covariance matrix of datapoints errors 
should be diagonal. Here we have used the same set of observational data, as in papers [HE] with some changes. 
The amplitudes and positions of 1st and 2nd acoustic peaks and position of first dip have been taken from results 
of WMAP mission team 7 , position and height of third peak have been taken from last recompilation of results 
of BOOMERanG experiment [Sj. At large angular scales (lower multipoles range) we use all datapoints from 
COBE experiment and WMAP pQ, except dipole. 

The CMB data are complemented by large-scale structure data: mass function and spatial distribution of 
rich galaxy clusters, temperature function of X-ray clusters, Ly-a forests of absorption lines in distant quasars 
spectra, peculiar velocities of galaxies. These LSS data can establish the amplitude and shape of spatial power 
spectrum of matter density perturbations. Determination of parameters has a problem of the degeneracy in 
dependence of observational manifestations upon the parameters. This problem is partially eliminated when 
we add additional measurements to the observational set: determinations of expansion dynamics by angular 
distance to Supernovae of la type, independent determinations of Hubble constant, constraints on baryon content 
from Big Bang nucleosynthesis theory. 

In sum there are 41 values in list of observational datapoints along with la statistical errors of their mea- 
surements. We consider all measurements to be independent and take their probability distribution as normal 
distribution. These values and errors are described in details in cited papers. Thus if we compare these results 
with those from paper |S] , we obtain the answer to the question "How improvements in measurements of CMB 
anisotropies affect the determination of parameters of observable Universe?" . 

LIKELIHOOD ANALYSIS 

Lets us have N observational characteristics and we are searching best-fit values of n cosmological parameters. 
In other words, we have the parameter set: 

P = (£lcDM,&de,tt u ,N v ,n b ,h,A s ,n s ,A t ,nt,T c ) 

and the set of observations: 

U = (A pi , £ Pl , A P2 , £p 2 , Ap 3 , £ P3 , £d 1 , Ci low , LSS, N S, h, SN) 

To find best-fit values for 11 parameters we are used the Levenberg-Marquardt algorithm |2j for minimization 
of function 

where jjj - observational value for some j-th characteristic, yj - its theoretically predicted value, Aijj - statistical 
uncertainty for measured value. Since the number of neutrino species N v is a discrete value, so we were searching 
for values of 10 parameters at N v held fixed at 1, 2 3. As in previous papers, instead of dimensional values 
of amplitude of scalar mode power spectra A s , we have used dimensionless value Sh, defined as r.m.s. of the 
fluctuation of matter density on the scale of present horizon of particle, relate as A s = 27r 2 ^(3000Mpc//i) 3+ " s . 

To find errors for parameters values one needs to carry out exploration of likelihood function profile C oc 
e~2* in parametric space to determine confidential ranges (marginalizing). The estimations of confidential 
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ranges for parameters are mostly based on the direct likelihood function integration ,9: , a kind of cumbersome 
computations. In papers we have proposed "economical" methods for marginalization procedure to avoid 
the direct integration. Here we propose a combined approach: likelihood function profile for parameter Xk is 
built using the minimization of % 2 in subspace of n — 1 parameters 

C(x k ) = e -3lx 3 (^ fc ^)-x 2 mi n] j ( 4 ) 

where a^/ _ best-fit values of cosmological parameters (i = 1, 2, . . . 12, i =/= k), for which x 2 -function has the 
minimum at fixed value of parameter Xk- This kind of representation for C(xk) can be obtained from the general 
integral form. Numerical experiment proves that results from proposed function C(xk) for confidential levels 
estimation virtually coincide with values, obtained by integration. 

RESULTS 

For historical reasons the common practice is to give a description of relic gravitational waves by the ratio of 
amplitudes of tensor T and scalar S contributions to the quadrupole {I = 2) component of CMB power spectrum. 
However, probability distributions for A t and T/S differ, and C{A t ) function is closer to the Gaussian shape 
than C(T/S). So we are using C{A t ) for determination of upper bounds and later recalculate it to T/S. We 
define the upper bound A\" at confidential level 2a as value for which the square under curve C(A t ) consist 
95.4% of the total square under this curve from to oo. 
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Figure 2. Likelihood functions C(A t ) = exp[— ^x(-^t)] f° r various sets of observational data. 

In terms of our technique CMB and LSS data have not sufficient statistical weight to make possible simulta- 
neous determination of both amplitude and slope of tensor power spectrum. In other words, if we would leave 
both A t and n t as free parameters, then A t could take arbitrary large value under condition n t — > — oo. If 
we keep the lower bound for n t fixed, then upper bound for A t depends on it. This problem finds its natural 
explanation by limitations of indirect method used. We have no data to distinguish the amplitude and slope 
separately. Fortunately, the majority of inflation models relate the n t with the slope of scalar power spectrum 
n s . So we have also analyzed the likelihood functions for the same observational data and some generic inflation 
models: with flat spectrum of tensor mode (n t — 0), natural inflation with n t — n s — 1 and chaotic inflation 
with n t = 0.5(n s — 1). The upper 2a constraints for them are following: A\ a = 1.9 • 10~ 5 for the first model and 
A^ 7 = 1.0 • 10~ 5 for the rest. Corresponding values for them are T/S — 0.6 and 0.18. These three models are 
more interesting from the standpoint of manifestations of tensor mode in the data of observations, so further 
analysis will include merely them. As far as likelihood functions C{A t ) for them are very close it is quite enough 
to analyze one of them, namely we take the model with n t — 0.5(n s — 1). 

In Fig. |21 it is shown how observational data influence the half- width of likelihood function. As we can see, 
adding the observational data with constraints on cosmological parameters themself, like Hubble constant, Big 
Bang Nucleosynthesis, SNIa observations and data on large-scale structure, decreases the level of confidence for 
the models with high tensor amplitudes. At the confidential level of 2a (95.4%) this amplitude cannot exceed 
~ 20% of scalar mode amplitude. In our previous estimations (see [§]), based on the data of balloon experiment 
BOOMERanG, this limitation was almost four times larger (at la C.L. the ratio was estimated as T/S = 1.7 
for the model with free n t ). So we clearly see the achievements of WMAP with high precision and covering over 
all sky. 
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The results for the whole parameter set are summarized in the Tabled There the constraints are tabulated 
for the case of chaotic inflation. As we can see, the accordance of these best-fit values with previous determina- 
tions E] is quite satisfactory. The new precise measurements of CMB temperature anisotropics significantly 
tightened confidential ranges and lowered upper constraints for f2„, t c and T/S. Obtained results agree with 
determinations of other authors, that used the WMAP data (those use somewhat different definitions). 



Table 1. Best-fit values of parameters, lower and upper bounds at 2a (95.4%) confidential level. 
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CONCLUSIONS 

From the standpoint of statistics interpretation of available observational data on large-scale structure of the 
Universe and CMB does not require the presence of considerable amount of relic gravitational waves. According 
to the basic assumptions of early Universe physics these relic waves have inescapable nature and common origin 
with matter density fluctuations seeding the large-scale structure. Within the framework of inflation theory the 
manifestations of gravitational background can be quite prominent, as long as amplitude of relic GWs directly 
connected to energy scales of inflation, e.g. Grand Unification Theory energy scale. Thus, the upper bounds on 
amplitude of GWs appear to be of the utmost importance for finding constraints on the moment and energy 
scales of inflation allowing to discriminate among models of inflation. 

The advances in measurements of CMB in space experiment WMAP substantially lowered upper bound for 
amplitude of tensor mode of perturbations (i.e. relic gravitational waves) to the level oiT/S < 0.6 (95.4% C.L.) 
for models with free slope parameter n t . For models with flat power spectrum of gravitational waves (n t = 0) , 
or some close to that spectra (n t ~ 1 — n s ), this limit appears to be the even lower T/S < 0.18 (95.4%). 



[1] http://lambda.gsfc.nasa.gov/ 

[2] Grishchuk L.P. et al. Gravitational wave astronomy: in anticipation of first sources to be detected // Usp.Fiz.Nauk.- 
2001.-171.-P. 3-59. 

[3] Durrer R. The theory of CMB anisotropies // Journal of Physical Studies.-2001.-5, N 2.-P. 177-215. 
[4] http://www.cmbfast.org 

[5] Durrer R., Novosyadlyj B. Cosmological parameters from complementary observations of the Universe // Monthly 
Notic. Roy. Astron. Soc.-2001.-324.-P. 560-572. 

[6] Durrer R., Novosyadlyj B., Apunevych S. Acoustic Peaks and Dips in the Cosmic Microwave Background Power 
Spectrum: Observational Data and Cosmological Constraints // Astroph. J. -2003.-583, N l.-P. 33-48. 

[7] Page L. et al. First- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Interpretation of the TT 
and TE Angular Power Spectrum Peaks // Astroph. J. Suppl Ser.-2003.-148.-P. 233-241. 

[8] Ruhl J.E. et al. Improved Measurement of the Angular Power Spectrum of Temperature Anisotropy in the Cosmic 
Microwave Background from Two New Analyses of BOOMERANG Observations // Astrophysical J.-2003.-599- 
P. 786-805. 

[9] Press, W.H. et al..// Numerical recipes in FORTRAN-1992.-New York:Cambridge Univ.Press. 

[10] Spergel D.N. et al. First- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Determination of 
Cosmological Parameters //Astroph. J. Suppl.Ser.-2003.-148.-P. 175-194. 

[11] Novosyadlyj B., Apunevych S Cosmological constraints on the amplitude of relic gravitational waves // Journal of 
Phys. Studies.-2005. (in press) 



(i 



